Distribution of parametric conductance derivatives of a quantum dot 
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The conductance G of a quantum dot with single-mode ballistic point contacts depends sensitively 
on external parameters X, such as gate voltage and magnetic field. We calculate the joint distribu- 
tion of G and dG/dX by relating it to the distribution of the Wigner-Smith time-delay matrix of 
a chaotic system. The distribution of dG/dX has a singularity at zero and algebraic tails. While 
G and dG/dX are correlated, the ratio of dG/dX and G(l — G) is independent of G. Coulomb 
interactions change the distribution of dG/dX, by inducing a transition from the grand-canonical 
to the canonical ensemble. All these predictions can be tested in semiconductor microstructures or 
microwave cavities. 
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Parametric fluctuations in quantum systems with a 
chaotic classical dynamics are of fundamental importance 
for the characterization of mesoscopic systems. The fluc- 
tuating dependence of an energy level Ej (X) on an ex- 
ternal parameter X, such as the magnetic field, has re- 
ceived considerable attention 0] . A key role is played by 
the "level velocity" dEj/dX, describing the response to 
a small perturbation In open systems, the role of 

the level velocity is played by the "conductance velocity" 
dG/dX. Remarkably little is known about its distribu- 
tion. 

The interest in this problem was stimulated by experi- 
ments on semiconductor microstructures known as quan- 
tum dots, in which the electron motion is ballistic and 
chaotic A typical quantum dot is confined by gate 
electrodes, and connected to two electron reservoirs by 
ballistic point contacts, through which only a few modes 
can propagate at the Fermi level. The parametric depen- 
dence of the conductance has been measured by several 
groups [^-js). In the single-mode limit, parametric fluc- 
tuations are of the same order as the average, so that one 
needs the complete distribution of G and dG/dX to char- 
acterize the system. Knowing the average and variance 
is not sufficient. Analytical results are available for point 
contacts with a large number of modes |p|-|l5|| . In this pa- 
per, we present the complete distribution in the opposite 
limit of two single-mode point contacts and show that 
it differs strikingly from the multi-mode case considered 
previously. 

The main differences which we have found are the fol- 
lowing. We consider the joint distribution of the con- 
ductance G and the derivatives dG/dV, dG/dX with 
respect to the gate voltage V and an external parameter 
X (typically the magnetic field). If the point contacts 
contain a large number of modes, P{G, dG/dV, dG/dX) 



factorizes into three independent Gaussian distributions 
|p|-|l2|. In the single-mode case, in contrast, we find 
that this distribution does not factorize and decays alge- 
braically rather than exponentially. By integrating out 
G and one of the two derivatives, we obtain the conduc- 
tance velocity distributions P{dG / dV) and P{dG / dX) 
plotted in Fig. ^ Both distributions have a singular- 
ity at zero velocity, and algebraic tails. A remarkable 
prediction of our theory is that the correlations between 
G, on the one hand, and dG/dV and dG/dX, on the 
other hand, can be transformed away by the change of 
variables G — (2e^//i) sin^ 9, where 9 is the polar coordi- 
nate introduced in Rcf. [Q. The derivatives d9/dV and 
dO/dX are statistically independent of 9. There exists no 
change of variables that transforms away the correlations 
between dG/dV and dG/dX. 

Another new feature of the single-mode case concerns 
the effect of Coulomb interactions |l6yi7t] . In the simplest 
model, the strength of the Coulomb repulsion is measured 
by the ratio of the charging energy e'^/G (with C the ca- 
pacitance of the quantum dot) and the mean level spacing 
A. In the regime e^/C 3> A, where most experiments 
are done. Coulomb interactions suppress fiuctuations of 
the charge Q on the quantum dot as a function of V or 
X, at the expense of fluctuations in the electrical poten- 
tial U. Since the Fermi level fi in the quantum dot is 
pinned by the reservoirs, the kinetic energy E = fj, — U 
at the Fermi level fluctuates as well. Fluctuations of E 
can not be ignored, because the conductance is deter- 
mined by E, an not by /i. An ensemble of quantum dots 
with fixed Q and fluctuating E behaves effectively as a 
canonical ensemble — rather than a grand-canonical en- 
semble. In the opposite regime e^/C ^ A, the energy 
E does not fluctuate on the scale of the level spacing. 
The ensemble is now truly grand-canonical. Fluctua- 
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tions of E on the scale of A can be neglected in the 
multi-mode case, so that the distinction between canon- 
ical and grand-canonical averages is irrelevant. In the 
single-mode case the distinction becomes important. We 
will see that the distribution of the conductance velocities 
is different in the two ensembles. (The distribution of the 
conductance itself is the same.) The difference between 
grand-canonical and canonical averages has been studied 
extensively in connection with the problem of the persis- 
tent current ]l8[-p0[ , which is a thermodynamic property. 
Here we find a difference in the case of a transport prop- 
erty, which is more unusual [ pTp^ . 

To derive these results, we combine a scattering formal- 
ism with random-matrix theory . The 2x2 scattering 
matrix S determines the conductance 

G = |5i2p, (1) 
and the (unscreened) compressibilities [T7| 

dE~2m dE' dX~2ni dX- 

(We measure G in units of 2e^ /h and Q in units of e.) 
Grand-canonical averages (• • ■)gc and canonical averages 
(• • ■)c are related by 
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The factor dQ /dE is the Jacobian to go from an average 
over Q in the canonical ensemble to an average over E 
in the grand-canonical ensemble. Conductance velocities 
in the two ensembles are related by 



(4) 



where \q and \e indicate, respectively, derivatives at con- 
stant Q (canonical) and constant E (grand-canonical). 
Derivatives dG/dV with respect to the gate voltage are 
proportional to dG / dQ in the canonical ensemble and to 
dG/dE in the grand-canonical ensemble. (The propor- 
tionality coefficients contain elements of the capacitance 
matrix of the quantum dot plus gates.) The two deriva- 
tives are related by 



dG dG rdQ 
dQ ~ dE \dE 



(5) 



The problem that we face is the calculation of the joint 
distribution of S, dS/dE, and dS/dX. In view of the 
relations (|)-(|) it is sufhcicnt to consider the grand- 
canonical ensemble. This problem is closely related to 
the old problem of the distribution of the Wigner- 
Smith delay times ti, . . . ,TAr, which are the eigenvalues 
of the N X N matrix —iS^dS/dE. (The eigenvalues are 
real positive numbers.) Interest in this problem has re- 
vived in connection with chaotic scattering p5|-p8[ | . The 
rates 7„ = 1 /r„ are distributed according to p8[ 
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FIG. y. Distributions of conductance velocities in a chaotic 
cavity with two single-mode point contacts [inset in (a)] , com- 
puted from Eq. ([id). Dashed curves are for /3 = 1 (time- 
reversal symmetry), solid curves for /3 = 2 (no time-reversal 
symmetry). (The case /3 = 4, which is similar to /3 = 2, is 
omitted for clarity.) The distribution of A9G/(?i? (grand- 
canonical ensemble) is shown in (a) and the distribution of 
dG/dQ (canonical ensemble) is shown in (b). (The conduc- 
tance G is measured in units of 2e^//i, the charge Q in units 
of e.) In (c) the distribution of XodG/dX is shown for the 
grand-canonical ensemble (the canonical case being nearly 
identical on a linear scale). The inset compares the canon- 
ical (C) and grand-canonical (GC) results for /3 = 2 on a 
logarithmic scale. 
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This distribution is known in random-matrix theory as 
the Laguerre ensemble, because the correlation functions 
can be written as series over (generalized) Laguerre poly- 
nomials psl. For = 1 we recover the result of Refs. 
p5| and ||27| . In our case N ~2. 

To compute the conductance velocities it is not suffi- 
cient to know the delay times t„, but we also need to 
know the distribution of the eigenvectors of the time- 
delay matrix —iS^dS/dE. Furthermore, we need the dis- 
tribution of —iS'^dS/dX. The general result containing 
this information is ESI 



P{S,TE,Tx) oc exp 



„ '1 vl 



X) 



X (detT^) 



-2^7V+3(/3-2)/2 



(7) 

dE" ' ^" dx" ■ 

The matrix te has the same eigenvalues as the time-delay 
matrix, but it is more convenient because it is uncorre- 
lated with S*, while the time-delay matrix is not. By 
integrating out te and tx from Eq. (|^) , we obtain a uni- 
form distribution for 5, as expected for a chaotic cavity 
pot - The resulting distribution of the conductance 
P{G) cx G~-'^+^/^, is the same in the canonical and grand- 
canonical ensembles, because S and dQ/dE are uncorre- 
lated [cf. Eq. (||)]. By integrating out S, tx, and the 
eigenvectors of te, we obtain the distribution (^) of the 
delay times. The distribution of tx at fixed te is a Gaus- 
sian. The scale of this Gaussian is set by the parameter 
Xq, which has no universal value [p2[. 



We are now ready to compute the distribution of the 
conductance velocities. Derivatives with respect to E 
and Q are related to the delay times by 



dG 
dE 
dG 

dQ 



c(ti - T2)y/G{l - G), 
Tl - T2 



2ttc- 



Tl + T2 
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(9a) 
(9b) 



where c G [— 1 , 1] is a number that depends on the phases 
of the matrix elements of S and on the eigenvectors of 
Te- Its distribution P{c) oc (1 — c^)~^+'^/^ is independent 
of Tl, T2, and G. The derivative dG/dX has a Gaussian 
distribution at a given value of S and te , with zero mean 
and with variance 



1 fdG 




= aTlT2 



where we have abbreviated a — 4A^/7r^Xg/3. Because 
the variance of dG/dX depends on dG/dE or dG/dQ, 
these conductance velocities are correlated. 

From the distribution (|^) of ri, T2, and the indepen- 
dent distributions of G and c, we calculate the joint 
distribution of G and its (dimensionless) derivatives 
Gx = XodG/dX, Ge = iA/2TT)dG/dE, and Gq = 
{\/2n)dG/dQ. The result in the grand-canonical and 
canonical ensembles is 



Pgc{G, Ge, Gx] 
PciG, Gq, Gx) = 
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{l-y)iP+3)/2^G(l-G)g{x) 



■ exp 



2(3 



exp 



G'x 



2(3x 



fix) = 8f3-'[xG{l -G) + 2G|], gix) = 8ix^l3)~'[G{l - G) - GU, Z = 3r'^-ir(/?/2)r(/3)r(3/3/2). 
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By integrating out G and one of the two derivatives 
from Eq. (p^, we obtain the conductance velocity distri- 
butions of Fig. |l|. (The case (3 — A\s close to /? = 2 and is 
omitted from the plot for clarity.) The distributions have 
a singularity at zero derivative: A cusp for [3 = 2 and 4, 
and a logarithmic divergence for (3 = 1. The tails of the 
distributions of dG/dX are algebraic in both ensembles, 
but with a different exponent, 

PGc{dG/dX) cx [dG/dX)-^--'^, (11a) 
Pc{dG/dX) cx {dG/dX)-^P-\ (lib) 

The distribution of DG/dE (grand-canonical ensemble) 



also has an algebraic tail [cx [dG / dE)^!^^"^], while the 
distribution of dG/dQ (canonical ensemble) is identically 
zero for \dG/dQ\ > tt. In both ensembles, the second mo- 
ment of the conductance velocities is finite for (3 = 2 and 
4, but infinite for (3=1. 

In conclusion, we have calculated the joint distribution 
of the conductance G and its parametric derivatives for 
a chaotic cavity, coupled to electron reservoirs by two 
single-mode ballistic point contacts. The distribution is 
fundamentally different from the multi-mode case, be- 
ing highly non-Gaussian and with correlated derivatives. 
(Correlations between G and the parametric derivatives 
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can be transformed away by a change of variables.) We 
account for Coulomb interactions by using a canonical 
ensemble instead of a grand-canonical ensemble. Our 
results for the canonical ensemble are relevant for the 
analysis of recent experiments on chaotic quantum dots, 
where the conductance G is measured as a function of 
both the magnetic field and the shape of the quantum 
dot 1^ . The grand-canonical results are relevant for ex- 
periments on microwave cavities p3| , ^ . Together with 
the theory provided here, such experiments can yield in- 
formation on the distribution of delay times in chaotic 
scattering that can not be obtained by other means. 

This problem was suggested to us by C. M. Marcus. 
We acknowledge support by the Dutch Science Founda- 
tion NWO /FOM and by the European Community (pro- 
gram for the Training and Mobility of Researchers). 
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